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1 Introduction 



While working at Los Alamos in 1943-44, von Neumann became convinced 
that the calculation of the flows of compressible fluids containing strong 
shocks could be accomplished only by numerical methods. He conceived 
the idea of capturing shocks, i.e., of ignoring the presence of a discontinu- 
ity. Employing a Lagrangian description of compressible flow, setting heat 
conduction and viscosity equal to zero, von Neumann replaced space and 
time derivatives by symmetric difference quotients. Calculations using this 
scheme were carried out; the approximation resulting from these calcula- 
tions (see ^0|) showed oscillations on the mesh scale behind the shock. Von 
Neumann boldly conjectured that the oscillations in velocity represent the 
heat energy created by the irreversible action of the shock, and that as Ax 
and At tend to zero, the approximate solutions tend in the weak sense to 
the discontinuous solution of the equations of compressible flow. 



In [14] it was counterconjectured that von Neumann was wrong in his 



surmise, i.e., that although the approximate solutions constructed by his 
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method do converge weakly, the weak limit fails to satisfy the law of con- 
servation of energy. 

In H J.Goodman and P.Lax investigated von Neumann's algorithm ap- 
plied to the scalar equation 

u t + uu x = (1) 

(it is called the Hopf equation), in the semidiscrete case. Using numerical 
experimentation and analytical techniques the demonstrated the weak con- 
vergence of the oscillatory approximations, and that the weak limit fails to 
satisfy the scalar equation in question. 

Von Neumann's dream of capturing shocks was realized in his joint work 



with Richtmyer in 1950, see [21|. Oscillations were eliminated by the judi- 
cious use of artificial viscosity; solutions constructed by this method converge 
uniformly except in a neighborhood of shocks, where they remain bounded 
and are spread out over a few mesh intervals. The limits appear to satisfy 
the conservation laws of compressible flow. The conservation of mass and 
momentum is the consequence of having approximated these equations by 
difference equations in conservation form; but the von Neumann- Richtmyer 
difference approximation to the energy equation is not in conservation form. 

In the paper [[n]] T.Hou and P.Lax compared the results of a von Neu- 
mann- Richtmyer calculation with the weak limit of of calculations performed 



by von Neumann's original method. P.Lax in the paper [1C] asserts that the 
difference scheme of von Neumann, because of the centering of the difference 
quotients, is dispersive; it is this quality that is responsible for the oscillatory 
nature of the solutions. 

E.Hopf H studied (|l|) defined the generalized solution for the this equa- 
tion. He considered the approximating equation 

ut + uu x = fm xx where fj, — ► (2) 

for the equation (U). 

By a generalized solution u of (||) or (|2|), [i = E.Hopf meant a function 
u that is measurable and quadratically integrable in every closed rectangle 
in the open semiplane t > and that satisfies the relation 

dxdt = 0, (3) 

where g is an arbitrary function of class C 1 in t > that vanishes outside 
some circle lying entirely in t > 0. He assert: Every limit function u obtained 
from the solution of (|2|) as fi — ► +0 is a generalized solution of ([!]). By a 
generalized solution u of (|2|) E.Hopf meant a function u that is measurable 
and quadratically integrable in every closed rectangle in the open semiplane 
t > and that satisfies the relation 

u 2 1 

uft + y/x + fJ-uf xx dxdt = 0, (4) 





u 2 1 


// 


u 9t + y3x 
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where / is an arbitrary function of class C 2 in t > that vanishes outside 
some circle lying entirely in t > 0. 

This method is called the "disappearing viscosity" method. It was de- 
veloped by E.Hopf U, O.A.Oleinik [||], P.Lax H], @. 

There is also the "zero dispersion limit" method developed by P.Lax 
|| , §§ Maslov V.P and his collaborators g, @. The idea is to use the 



following equation 

u t + uu x = e 2 u xxx , e -> 0. (5) 

for the (|l|). V.P.Maslov and his collaborators constructed an asymptotic 
solutions for the (|l|) . They speculated on the fact that equations @ and (||) 
have solutions in the class C°° functions. For example, a particular solution 
of the equation (||) is the function 

9 fx — ct\ 
c + 4 - 12 cosh -2 , c > 0, 



which represent so-called an infinitely narrow soliton. A particular solution 
of the (p[) will be, for instance, 



2 tanh 



x-2t 



which converges to a discontinues function. 

In the other hand, in the paper |30[ by S.L.Sobolev (1936) were in- 
troduced mathematical basics of the theory of generalized functions, was 
developed the idea of a generalized function as a functional and was intro- 
duced the concept of generalized solutions of a linear differential equation. 
These generalized functions and generalized solutions were developed by 



L.Schwartz [29|. However, Sobolev-Schwartz distributions can not applied 
to nonlinear differential equations. For example, to substitute generalized 
function H(x — vi) (where H is the Heaviside function) into the ([!]) one 
need to define the product of two distributions H and H' for the term uu x . 
However, in 1954 L.Schwartz showed that it is impossible to introduce an 
associate multiplication in the space of distributions. 

Starting from 1982 in the works by J.-F. Colombeau ||], M. Oberguggen- 
berger p|], H. Biagioni [§, E. Rozinger f2l, A.Y. Le Roux jl, Yu. Egorov 



[§, J.-A. Marti 0, A. Delcroix, D.Scarpalezos §, B. Keyfitz 0, A. An- 
tonevich, Ya. Radyno |l|], T. Todorov |3l|], S. Pilipovic |1| and others, a 
new theory of generalized functions is developed. Such functions form the 
algebra and contain distributions. 

In general, nonlinear generalized functions are classes of equivalent 
smooth functions. Clearly that, one should pay attention to this approach 
in order to consider nonlinear differential equations. However we are now in 
a position to develop a new point of view on generalized functions and their 
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applications to nonlinear equations. Namely, it is necessary to use an integral 
nature of a conservation law. Conservation laws are integral expressions 
from physical point of view and it is natural to consider an integral form 
of conservation laws. Moreover, we want to develop new point of view 
on conservation laws using the concept of functionals with values in the 
Non- Archimedean field of Laurent series. We call such functionals as Re- 
distributions [25], [26]. In addition, we give the definition of the special kind 
of solutions of the some conservation laws in the sense of R(e)-distributions 
and consider the method for the numerical calculations of the smooth shocks 
and soliton like solutions of the Hopf equation and equations of elasticity 
theory in the mentioned sense. This method based on orthogonal system of 
the Hermite functions as a base for calculation of such solutions (i.e. shocks 
and infinitely narrow solitons). Calculations of profiles of infinitely narrow 
soliton and shock wave are reduced to the nonlinear system of algebraic 
equations in R n+1 , n > 1. We proved, using the Schauder fixed point 
theorem [28], that the mentioned system has at least one solution in R n+1 . 



We showed that there is possibility to find out some of the solutions of this 
system using the Newton iteration method [Hfl . We considered examples 
and numerical tests. We also should emphasis that proposed numerical 
approach do not use a difference scheme. 

First, let us consider a bit of theory which we will apply to conservation 
laws. 



2 Non- Archimedean field of Laurent series and 
R(e}— distributions. 

The theory of Non- Archimedean fields was considered in the book by 
A.H.Lightstone and A.Robinson flTq]. 

Definition 2.1 A Laurent series is a formal object 



n=0 



where k is a fixed (i.e., fixed for this Laurent series), each £j 6 R, and either 
£fc ^ or each £j = 0. 

The Laurent series ^JS=o £n£ n j where £o = 1 and £ n = if n > 0, is 
denoted by 1. It is easy to see that the Laurent series is a field. Let us 
denote it by R(e). The norm on the field of Laurent series can define 

\x\ u = e~ u ^ for each x G R(e) 
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(in place of e can use any number greater then 1). The function v(x) is a 
Non-Archimedean valuation. Define 

(oo \ oo 

J2 tn+kE n+k j = k if J2 Zn +k e n+k + 0, 6/ 0. 
n=0 / n=0 

The norm | • \ v have properties 

1. = if and only if x = 0, 

2. |xy|j, — \^\u ' \y\i/i 

3. \x + y\ u < max{\x\ u , \y\ u } . 

Here, we propose a general construction of the R(e)-valued generalized func- 
tions [24]. These objects are a natural generalization of Sobolev-Schwartz 
distributions. We call them as R(e)-distributions. 

1. Consider all functions f(x,e) 6 C°°(R x (0, 1)) such that integrals 

+oo 

j f(x,e)ip(x)dx 

— oo 

exist for any e and for all ip{x) from a given class of functions X (X 
can be Cg°(R), S(R) and etc.). 



2. Suppose also that J f(x, e)^(x)dx is a number a/ j£ (^) from the field 

— oo 

of Laurent series R(e). 

3. The two functions f(x,e) and g(x,e) call equivalent with respect to 
test functions X if and only if 

+oo +oo 

J f(x,e)ip(x)dx = a /)£ (V>) = a gt£ (ip) = J g(x 1 e)^(x)dx. 

— oo — oo 

The equality means in sense of the field of Laurent series R(e) for all 
functions ip G X. Classes of equivalent functions call R(e) -functions. 
The expression 

+oo 

J f(x,e)ip(x)dx 

— oo 

associates a number from R(e) with every if). Such a quantity is called 
a functional. In this linear functional map X into the Non- 

Archimedean field R(e). Call these functionals as R(e) -distributions. 
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Thus, 

Proposition 2.1 R(e) -function f(x,e) = if and only if 

+00 



J f(x,e)ip(x)dx = E R(e) 



for every ij) from X. 

The set of all R(e)-distributions denote by %{%) 

Remark 2.1 Recall that the idea of representation of a function f 6 
Lf (R) in terms of a linear functional 



i-oo 

C °°(R)9^ j f(x)^(x)d 



, 'x € R 

—00 

+00 

feri „„ ™«-*„ ™ WS!(t( ,n *tf „ / 6 LL(R) and / / (a # W(i x = 

— OO 

for any ip £ Cq°(R) then f = almost everywhere. 

Let us consider an example of the R(e)-distribution. 

Example 2.2 Take X = C§°(R) and f(x,e) = <p{x/e), (p(x) £ Cg°(R) 
then R(e) -distribution can write in the following form. 

+00 +00 +00 

V>'(o) 



ip(x/e)ip(x)dx = e J (p(x)dxt()(0) + e J xip(x)dx — — h . . . . 

—00 —00 
Note that ip(x/e) converges to the function 

UUC ,-,X°). tf* = o, 



0, ifx^O. 

Last function almost everywhere equals to zero. 

Like Sobolev-Schwartz distributions we can differentiate R(e) -distribu- 
tions. For example, 

+00 +00 

—ip(x/e)ip(x)dx = — cp(x/e)— — ip(x)dx, 

dx / dx 



(p(x/e)-^-ip(x)dx = — e [ (p(x)dxtf)' (0) — e 2 [ xip(x)dx^ ' ' 



dx^ y ' J ' \ / j t \ j 1 , 

—00 —00 

It is evident that R(e)-distributions are more general objects than Sobo- 



lev-Schwartz distributions | 29fl , [ 30fl . 
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3 Conservation laws. Non- Archimedean approach. 

A conservation law asserts that the rate of change of the total amount of 
substance contained in a fixed domain G is equal to the flux of that substance 
across the boundary of G. Denoting the density of that substance by u, and 
the flux by /, the conservation law is 

±fu(t, = /.«& 
dt J G JdG 

Applying the divergence theorem and taking d/dt under the integral sign we 
obtain 

/ (u t + div f)dx = 0. 
JG 

Dividing by vol (G) and shrinking G to a point where all partial derivatives 
of u and / are continuous we obtain the differential conservation law 



ut(t, x) + div f(u(t, x)) = 0. 

Note, that if f(u) = u 2 /2 then we obtained the Hopf equation ([l]). In 
general, previous calculations based on the following well known proposition. 

+oo 

Proposition 3.1 If G G L[ 0C (R) and J G(x)ip(x)dx = for any ip G 

— oo 

Cg°(R) then G = almost everywhere. 

Definition 3.1 Let us consider two sets of the smooth functions, depending 
on a small parameter e G (0, 1]. Let us take all functions v(t,x,e) which have 
the type 

fx — ct 

v(t, x, e) = Iq + Alip 



e 

+oo 



Iq, A/, c are real numbers, Al ^ and ip G S(R), J f{y)dy = 1. We denote 

— oo 

this set of functions by I. We call I as a set of infinetely narrow solitons. 

Definition 3.2 Now, let us take all functions w(t,x,e) which have the type 

x — at 

w(t, x, e) = ho + AhH ' 



x -|-oo 

ho, Ah, a are real numbers, Ah ^ and H(x) = J 6{y)dy, j 9{y)dy = 1 

— oo — oo 

and 6 G §(R). We denote this set of functions by J. We call J as a set of 
shock waves. 
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It is natural to consider conservation laws as an integral expressions 
which contain the time t as parameter. Therefore, we introduce the following 
concept. 

Definition 3.3 The function v € I (or w £ J) will be a solution of the 
Hopf equation up to e~ p , p £ No in the sense of R(e) -distributions if for 
any t <G [0, T] 



it +00 
/ {v t (t,x,e) +v(t,x,e)v x (t,x,e)}ip(x)dx = ^2^ k e k £K(e), (6) 
-00 k =P 

+2° +00 
/ {wt(t,x,£)+w(t,x,£)w x (t,x,£)}ip(x)dx = '^2r] k £ k £'R(e) (7) 

-00 k =P 

for every -ip £ S(R). In case when p is equal to +00 the function v(t,x,£) 
(or w(t,x,£)) exactly satisfies the Hopf equation in the sense of Re- 
distributions. 

Certainly, one can consider instead of the Hopf equation some conserva- 
tion law. 

Prom mathematical point of view, we deal with a infinitely differentiable 



functions in definitions 3.1 and 3.2, so that we avoid the problem of distri- 
bution multiplication. From physical point of view, functions from the set 
I or J can describe fast processes. Mathematical models of such processes 
based on functions from I or J may give additional information and take in 
account a short zone where physical system make a jump from one position 
to another. 

Thus, we will consider solutions of the Hopf equation which are infinitely 
narrow solitons or shock waves. It easy to see that 

v(t,x,e)^\ jo + AMO), if x = ct, ag £ ^ Q 
i h, if ct. 

w(t,x,£) — > ho + AhH(x — at), as e — > 
H is Heaviside function. 

4 Method for the numerical calculations of the mi- 
croscopic profiles of soliton like solutions of the 
Hopf equation in the sense of R(e}— distributions. 

Thus, conservation laws are integral expressions. Therefore, it is natural, 



that one can interpret the Hopf equation in the sense of the definition 3.3 



S 



We will seek a solution of the Hopf equation in the type of infinitely 
narrow soliton, i.e. let us v £ I. Substitute v(t,x,e) into integral expression 
(0) using the following formulas 



+00 
./ <)t 



-00 



I L f^\\^ x )dx = 2^ fe+1 m4v (fc+1) (ci) ) (8) 

^ ^ £ ' ' k=0 



8 ( 1 2 / X — ct 



dx 1 2^ 



) U(x)dx = Y,- £k+1 9k^ (k+1) (ct). (9) 

' ' fc=0 



We denote 



r X 

m fc (^)= / y«<p{y)dy, 5fc (<p) = / y k ^-dy, k = 0,1,2,... . (10) 



Thus, we obtain 

+00 



/ {v t + vv x }^dx = ^{Al(c-l )m k -(Al) 2 g k }e k+1 ^ 1 -^. (11) 



V>( fc+1 )(ci) 

fc=0 

From the last expression we have conditions for the function <p(x). Namely, 



9k(<p) - ?-£f-Tnk(<p) =0, A; = 0,1,2... . (12) 



From the first (k = 0) we have 



-oc 



c ~ l ° = 1 / ^(x)^. (13) 



A/ m 2 

—00 

Hence, we can rewrite conditions (|l2|) as follows. 

+00 +00 +00 

tp 2 (x)dx- / x k ip(x)dx= / x k ip 2 (x)dx, k = 0, 1,2... . (14) 



Now, let us prove the following lemma. 

Lemma 4.1 For any non-negative integer n exists such function ip E §(R), 
99 ^ which satisfies the following system of non-linear equations: 

+ OO +OO +CX! 

„k , 2 / 



x K ip(x)dx = / x K ip' s (x)dx/ / ip z {x)dx fe = 0,l,2...n. (15) 
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Proof. First, we will seek function (p(x) in the following type: 

(p{x) = coh (x) + cihi(x) + . . . + c n h n (x), (16) 

where 

hb,(x) = , e~ x I 2 are Hermit functions. (17) 

Then we substitute the expression ( |l6l) into conditions (^). After that 
we will have nonlinear system of n + 1 equations with n + 1 unknowns 
(co, c l) c 2i • • • ) c n ). We write this system by the following way. 



Ax = 'N(x), x = (co,ci, . . . ,Cn) (18) 
A is a matrix with elements 



-l-oo 

A kj = J x k h j {x)dx = (-i) j i k V2^hf\o), i = k,j = 0, 1, 2, 



n 



N is nonlinear map such that 

^ n n 

= — p 5^(JV(fc)^ f)eit = ^ / fc (x)e fc (19) 

" X " k=Q k=0 

Vector e*fc = (eo, &\, . . . , e n ) such that = 1 and = for all j / k. N(k) 
are matrices with elements 

+oo 

Nij{k)= j x k h i (x)h j (x)dx, i, j,k = 0,1,2 ... n (20) 

— oo 

and functions 



x 



Note that functions fk(x) are continuous everywhere except x = and 
|/fcO*OI < 11-^(^)11 due to Cauchy-Bunyakovskii unequality. Matrix A is in- 
vertible for any n because of det(A) is a Wronskian for the linear independent 
system of Hermit functions ho(x), hi(x), ... h n (x) and 

(n+l) 

det(A) = (2ir) — W(h (0), h^O), . . . h n (0)). 



We can write the system ( |I8| ) as 

n 

fk@)A-% = F(x) or x = A-^Nix )) = F(x). (21) 



x = 

k=0 
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Let us describe the function F : R n+1 i — ► R n+1 . It is continuous except 
x = and bounded. Indeed, 



n 



\\F(x)\\ < WA^W £ \\N(k)\\, r n = WA-'W £ W(k)\\. (22) 

fc=0 fc=0 

Let us consider function J4(x). It is continuous function everywhere in R™ +1 
except x = and, moreover, N(R n+1 \{0}) C LTi where EEi = {z € R ra+1 : 
zq = 1} is a plane. Further 74 _1 (Ili) = LT2 where II2 = {y G R n+1 : 

n 

^2 a 0jUj = 1} is another plane. 

fc=0 

+00 

A 0j = J hj(x)dx = (-i) j V2^hj(0), i = v 73 !, j = 0,1,2,... n 



-00 



Thus, we can consider the function F(x) which is defined on the convex 
compact set C n = n 2 P| -^[0, ?~n] such that F : C n 1 — > C n , where B[0,r n ] is 
a closed ball with radius r n . Function F is continuous on the C n because of 
C n . Now we can use J.Schauder theorem. 

Theorem 4.2 (Schauder fixed-point theorem p8|] ) Let C be a com- 
pact convex subset of a normed space E. Then each continuous map F : 
C 1 — ► C has at least one fixed point. 

Hence, we can conclude that our system (^Tj) and therefore system ( |T8| ) 
has at least one solution. Thus, there is a function <p(x) which satisfy to 
conditions (FL5|) proposed lemma. 



Remark 4.1 Let us a function <p(x) satisfies lemma condition. If (5 G R 
then the function (p(x + (3) also satisfies lemma condition. Moreover, if 



00 

2 



if (x) dx = a then if (ax) satisfies lemma condition. 



Thus, we can formulate the following result. 



Theorem 4.3 For any integer p there is a infinitely narrow soliton type 
solution of the Hopf equation (in the sense of the definition \3.3{ ) up to e~ p 
with respect to the norm \ ■ \ v , i.e. 



( x — ct\ 



v(t,x,e) = l + Al<p[ , (23) 
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-TOO 

Iq,AI,c are real numbers, Al ^ and ip £ S(R), J <p(y)dy = 1. Moreover, 



c-k 
Al 



(p 2 (x)dx. 



(24) 



For example, calculations in case p = 7 give the "profile" <p(x) (see Fig. 
|]) for the infinitely narrow soliton v(t,x,e) = <p {j-^)'- 



<p{x) - j C ° | c 2( 4a;2 ~ 2 ) + c 4 (16x 4 - 48x 2 + 12) | e _ a: 2 /2 



(25) 



where c = 0.66583, c 2 = -0.23404, c 4 = 0.05028, c = 0.25032 (c is a 
velocity of the soliton). Numbers Co, C2, c 4 and c were found approximately 
by iteration method using the following sequence. 

x m+1 = A-\JHx m )) t m = 0, 1, 2, . . . . (26) 

Matrix A and a nonlinear !N were introduced in the lemma proof. 




0.8 
0.6 
0.4 
0.2 

-0.2 



Figure 1: The case p = 7, c = 0.25032. Figure 2: The case p = 13, c = 0.35442. 

Calculations of soliton-like profiles <p(x) for the Hopf equation in case 
p = 13, 15, 17, 19, 21 give us pictures (Fig. g |,|, |, |). 

For the p greater than 21 matrix A is close to singular and calculations 
can be inaccurate. 



Calculations of the microscopic profiles of the 
shock wave solutions of the Hopf equation in 
the sense of R(e)— distributions. 



A solution of the Hopf equation in this case we will seek in the set J. Namely, 



12 




(X — Clt 
— - — 

ho, Ah, a are real numbers, Ah ^ and 

x +00 

K(x) = J e{y)dy, J e(y)dy = 1, 6 G §(R). 

—00 —00 

Substitute w(t, x,e) into the integral expression @ using the following 
formulas 



7s { A ' (^) } * {x)dx = D-^t- (27) 



-00 
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We denote by 

+00 

m k {9)= [ y k 6(y)dy, r k (9) = I x k 9(x) I / 0{y)dy I dx, £ = 0,1,2, 



(29) 



Thus, we get 

/ {w t + ww x }ipdx = J2{( Ah ) 2r k- Ah(a-h )m k }e k ^ j ; . (30) 

fc=0 



From the last expression we have conditions for the function 9{x) 

r k (0) - ^^m k (e) = Q, k = 0,1,2... . (31) 



From the first (A; = 0) we have 

o 

9{x) I / 9(y)dy I dx = (32) 



+00 / z 

a — ho f „, s / /" „ , v , \ , I 
A/i = 

Therefore, we can rewrite ([H]) as 
+00 +00 
i y x fe 0(x)dx= y a; fc 0(x) I / (//),/// I r/.r /,- 0.1.2... . (:«) 

—00 —00 

The same method one can prove that there is such function 9(x) £ §(R) 
which satisfies the following conditions 

+00 +00 / X \ 

1 1 x k 9(x)dx = y x k 9(x) I y 9{y)dy J dx k = 0, l,2...n. (34) 
—00 \— 00 / 

Thus, we can formulate next result. 

Theorem 5.1 For any integer p there is a shock wave type solution of the 
Hopf equation (in the sense of the definition S.c) up to e~ p with respect to 
the norm \ ■ \ u . 

(x — at \ 
£ j ' ^ 

x +00 

/io, A/t, a are real numbers, Ah 7^ and K(x) = j 9(y)dy, j 9(y)dy = 1 

—00 —00 

and 9 £ §(R). Moreover, 

a -h 1 

~A/T = 2' (36) 



14 



Note that the condition fl3q ) is Rankine — Hugoniot condition for the 
velocity of a shock wave. 

As in previous section we seek function 9{x) in the following type: 



(p(x) = a ho(x) + a\hi(x) + . . . + a n h n (x), 



(37) 



where h^{x) are Hermite functions. Calculations in case p = 7 give the fol- 
lowing "profile" 
h = 0,Ah= 1) 



x — at 

lowing "profile" {K(x)) for the shock wave w(t,x,e) = K ( ) (where 



K(x) 



X , 



.„ + c 2 (4^- 2) + o(16 T - -4 8T ^ + 12) e . T , /2<jT 



where c = 0.79617, c 2 = -0.53004, c 4 = 0.17923, c = 1/2 is a velocity of 
the shock wave (see Fig. ^). Numbers cq, C2, C4 were found approximately. 

Note that the function K(x) is not unique. There is a different function 
K\(x) which satisfies mentioned above conditions. It has the following type 



K, (x) 



+ 



cp ci2r c 2 (4r 2 -2) 1 ^ /2 

v 7 ^ vW^f vW^f r 



c 3 (8r 3 - 12r) c 4 (16r 4 - 48r 2 + 12) 



(39) 



where c = 0.18357, ci = -0.73567, c 2 = 0.74733, c 3 = 0.15327 c 4 = 
—0.29539, c = 1/2 is a velocity of the shock wave (see Fig. ||). Coefficients 
Co, ci, c 2 , C3, c 4 were found approximately by the Newton iteration method. 
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Figure 7: Graph of the function K(x). Figure 8: Graph of the function Ki(x) 



Taking in account the Rankine — Hugoniot condition ( p6[) we also have 
graphs (Fig. ||, 10) as a shock profiles. 
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Figure 9: First shock profile 1 - K(x). Figure 10: Second shock profile 1 — Ki(x). 

Here we describe how it is possible to find coefficients cn, c±, . . . , c n in this 
case by the Newton iteration method for the following system of nonlinear 
equations. 

n 

P(c)=A&-2j2{S(k)c,c)e k =0, c=(c ,c 1 ,...,c n ) (40) 

k=0 

Vector e k = (eo, ei, . . . , e n ) such that e k = 1 and ej = for all j / k. S(k) 
are matrices with elements 

+oo X 

Sij(k) = J x k hi(x) ■ j hj(y)dydx, k = 0, 1, 2 . . . n (41) 

— oo — oo 

Matrix A have elements 

+oo 

Aij = / x t hj(x)dx, i, j = 0,1,2 ... n (42) 



-oo 



We can write the formula for the Newton iteration method [12]. 



Xm+l = X m ~ [P'(x m )] 1 [P(Xm)\ , (43) 

where [P'(x)] is a linear map depending on the vector x. 

[P'(x)] [h]=Ah-2l J2(S(k)x, h)e k + ^2(S T (k)x, h)e k \ (44) 

{k=0 k=0 J 

Calculations of shock profiles K(x) for the Hopf equation in case p = 8, 
9, 10, 11, 12, 13 give us the following pictures (Fig. [□], 0, ||, 0, ||, 0). 
Here, we show only two different types of the shock type solutions of the 
Hopf equation. We can find more solutions if we take a different initial data 
for the Newton iteration method. 
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Figure 11: Shock profiles when p = 8. Figure 12: Shock profiles when p = 9. 
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Figure 13: Shock profiles when p = 10. Figure 14: Shock profiles when p = 11 
Remark 5.1 It is not easy to see that there is exist function 

oo 

0{x) = ^2 a n h n (x), a = (a , oi, . . . , a n , . . . ) G h, 



n=l 



such that 



^ J x k 9{x)dx= J x k 6{x)[ j 6{y)dy J dx, k = 0,1,2... . 



(45) 



(46) 



— oo — oo 

We think that it is true. 



6 Calculations of the microscopic profiles of the 
shock wave solutions of equations of elasticity 
theory in the sense of R(e)-distributions. 



Let us consider the following system. 

ut + (u 2 ) x = a x (the conservation law for momentum) 
a t + ua x = k 2 u x (the Hooke law) 



(47) 
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Figure 15: Shock profiles when p = 12. Figure 16: Shock profiles when p = 13. 

Here, u is the velocity of a medium and a is the stress. We suppose that 
density of a medium is equal to the constant 1 and k 2 some constant. 
We will seek for a solution of this system in the following form 

fx — vt \ 

u(t, x, e) =u + AuU , (48) 



x +oo 

uq, An, v are real numbers, An ^ and U(x) = J U(y)dy, J U{y)dy 

— CO — CO 

and U G 8(R). 



x — vt 

a{t, x, e) = a + AaS ( ), (49) 



X -|-0O 

o"o, Act, v are real numbers, Act / and £(x) = J Y^(y)dy, J T,(y)dy = 1 

— oo — oo 

and £ € §(R). Note that v is a velocity of the shock waves. 
In the other hand, we suppose 

U(x) = aoh (x) + aihi(x) + . . . + a n h n (x), a = (a , ai, . . . , a n ), (50) 



E(x) =c h (x)+c 1 hi (») + ... + c n h n {x), c = (c , ci, . . . , c„) (51) 

where hk{x) are Hermite functions. 

We understand the solution of the system in sense of R(e)-distributions. 

Definition 6.1 Functions u G J and a £ J is a solution of the system ftffi ) 
up to e~ p , p £ Nq in the sense of R(e) -distributions if for any t £ [0, T] 
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•y +00 

/ {u t (t, x, e) + 2u(t, x, e)u x (t, x, e) - a x (t, x, e)} ip(x)dx = ^ £ k e k 6 R(e), 



k=p 

(52) 



+00 



J {a t (t, x, e) + u(t, x, e)a x (t, x, e) - k 2 u x (t, x, e)} ip(x)dx = ^ r] k e h G R(e) 
-00 fc =P 

(53) 

for every ijj £ §(R). 

In case when p is equal to +00 functions u(t, x, e) and <r(t, x, e)) exactly 
satisfies the system (47) in the sense of R(e)-distributions. 



Substituting u and a into (|52|), (|53| ) we get the following relations for 
the moments. 



{2u Au - vAu}m k (U) + 2(Au) 2 m fc ([/[/) - A<jm fc (£) = 0, fc = 0, 1, 2, . . . n 

(54) 

{uoAcr - wAa}m fc (S) + AuAam k (Y,U) - k 2 Aum k (U) = 0, k = 0, 1, 2, . . . n 

(55) 

We denote as usual by 

+00 

m k (U) = / x k U(x)dx, (56) 



OO / X 

m k (UU)= [ x k U(x) I / U(y)dy I dx, k = 0,1,2,... , (57) 



CXI 

m fc (Si7) = / x fc S(x) I / C/(y)dy I dx, k = 0, 1, 2, . . . (58) 



It is easy to find u from ( p4j ) when k=0. Indeed, 

{2u Au - vAu} + (Au) 2 - Act = 0. 
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Therefore, 



Act . , 

v = 2u + Au-—. (59) 
An 



Substitute v into the (Q). We have 

{Act - (Au) 2 }m k (U) + 2{Au) 2 m k (UU) - Acrm fc (S) = 0, ft = 0, 1, 2, ... n 

(60) 

Because of AU and Act some real numbers, therefore, all three vectors 
with coordinates m k (U), m k (UU) and m k (H), fc = 0, 1, 2, . . . n, respectively 
should be collinear. However, 

mo(U) = mo(S) = 1. 

Hence, 

mjt(Z7) = mjfc(S), = 0, 1,2, . . . n. 
Thus, 5 = c and from ( ]60j ) follows that 

m k (U) = 2m k (UU), k = 0, 1, 2, . . . n. 

This system we already know how to solve by the Newton iteration method. 
See conditions (p3) and solution in this case. 



Substitute v into the (55) and take into account previous equalities we 
have 

^ u Aa - AnAcrl m fc (S) + AuAam k (^) - k 2 Aum k (Y,) = 



An 



(61) 



where k = 0, 1,2, ... n. The last expression gives us relation for constants 
Act, Au, uq, k 2 . Namely, 

(A ' Ji uoAct — AnAcrl + -AuAa — k 2 Au = (62) 

or 



An 2 



(Act) 2 - ^n + ^An j AuAa - k 2 (Au ) 2 - (6:i) 
If An, no, k 2 are known then from the last equation one can find Act 



Act 1j2 = 1 f n + ^An j ± ~|An| J (n + ^An) 2 + 4A; 2 (64) 
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In particular, if An = — 1, uq = 1 then 



Aai : 



4 4 



16fc 2 . 



Shock profiles of the considered system (| 
(Fig. 



3 E 

k 2 = 



19, Eq 



one can find on pictures 
We considered case when p = 13 and An = — 1, 
no = l, le" = 0.1. We can also take any real o"o but here we took oq = 0.5 
and then calculated Aa\ (Fig. [L7], [L9|, the velocity of shocks is v = 1.1531 ) 
and A<72 (Fig. |(], the velocity of shocks is v = 0.34689). We consider 
two different types of shock profiles. The first is on the Fig. 0, [[§ The 



second is on the Fig. 19, 20 



Theorem 6.1 For any integer p there is a solution of the system of equa- 
tions fflj ) in the sense of the definition \6.1\ Moreover, 



2uq + An 



An 



and 



(Aa) 



n + ^An ) AuAa - k 2 (Au) 2 



0. 
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Figure 17: Shock profiles of velocity and Figure 18: Shock profiles of velocity and 
stress. stress. 

Let us consider the following system. 

pt + (pu) x = (the conservation law for mass) 

(pu)t + (pu 2 ) x = a x (the conservation law for momentum) (65) 
at + ua x = k 2 u x (the Hooke law) 

Here, n is the velocity of a medium and a is the stress. We suppose that k 2 
is some constant. 

Definition 6.2 Functions u £ J , p £ J and a £ J is a solution of the 
system ([^) up to e~ p , p G No in the sense of R(e) -distributions if for any 

te[0,T] 
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Figure 19: Shock profiles of velocity and Figure 20: Shock profiles of velocity and 

stress. stress. 



/ {p t {t,x,e) + p x u + pu x }i){x)dx = ^£ k s k e~R(e), (66) 



+S° +oo 

/ {p t u + pu t + p x u 2 + 2puu x -a x }^(x)dx = ^2( k e h eK(e), (67) 

-oo k =P 



/ {a t + ua x - k 2 u x } tp(x)dx = £ R(e). (68) 

-oo k =P 

for every tp £ S(R). 

In case when p is equal to +oo functions u(t, x, e), p(t, x, e) and a(t, x, e) 



exactly satisfies the system ffiq) in the sense of R(e) -distributions. 



We will seek for a solution of this system in the form fl4S|), (|49|), ((50), 



x — vt 

p(t,x,e)= Po + ApR[ ), (69) 



po,Ap,v are real numbers, Ap ^ and R(x) = j R(y)dy, J R(y)dy = 1 

— oo — oo 

and R £ §(R). We suppose 

R(x) = b h (x) + hhi(x) + . . . + b n h n (x), b = (6 , h, ■ ■ ■ , K), (70) 



22 



Substituting fl4S| ) and (^) into (q^) we get the following relations for 
the moments. 



{-vAp + Apu }m k (R) + ApAum k (RU) + p Aum k (U) + ApAum k (UR) = 

(71) 

where & = 0, 1, ... n.. Note that mo(RU) = 1 — mo(UR). Suppose k = 
and we get 

—vAp + Apuo + ApAu — pqAu = 

or 

Au 

v = uq + Au + pq— (72) 
Ap 

Last expression gives us the following 

-Au(p + Ap)m k {R) + ApAum k (RU) + poAum k {U) + ApAum k (UR) = 

(73) 

where /c = 0, 1, ... n.. 

All four vectors with coordinates m k (R), m k (RU), m k (U), m k (UR) 
should be collinear. Consider m k (R) and m k (U). Because of mo(R) = 
m o(U) = 1 then m k (R) = m k (U) for k from to n and therefore a = b. 
From last equality we have 

m k (RU) = m k {UR) = m k {UU) 
where k = 0, 1, . . . n.. Thus 

-Au(p + Ap)m k (U) + 2ApAum k (UU) + p Aum k (U) = (74) 

where k = 0, 1, . . . n.. It means 

m k (U) = 2m k (UU), fc = 0, 1, . . . n. (75) 



Substituting (|48|), (^9|) and (|69|) into (|67|) we get the following relations 
for the moments. 



u Ap(u - v)m k (R) + p Au(2u - v)m k (U)+ 

+AuAp{2u - v){m k (UR) + m k {RU)}+ . . 

+2p (Au) 2 m k (UU) + p{Au) 2 m k (RU 2 )+ K ' 

+2Ap(Au) 2 m k (URU) - Acrm fc (S) = 

where k = 0, 1, . . . n.. Suppose k = and we get 

Po(Au) 2 Q^ + l)+Aa = 0. (77) 
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Moreover m&(i?) = mfc(S), k = 0, 1, . . . n and then b = c. Finally, 



b = c and mk{R) = 3mfc(C/[/ 2 ), k = 0, 1, 



Substituting (g|) and (|g) into (| 
the moments. 



we get the following relations for 



{u Acj - t;Acr}m fc (S) + AuAam k (EU) - k 2 Aum k (U) = (78) 
where k = 0, 1, ... n.. Suppose k = and using expression for the velocity 



(ED we get 



AaAu (j^ + + A; 2 An = 0. 

From the (|7^) we can also find the following equality for the velocity 

1 , 2 An 
v = n + Au + -n - fc ^ 

It is well known result in the elasticity theory. 

Thus, if pq, Au and k 2 are known then the rest constants we can find 
from the system 



(79) 



(80) 



Po(Au) 2 (p + Ap) + AaAp = 0, 
Act (po + 3A/9) + /c 2 Ap = 0. 



Hence, 



and 



Act 



2k 2 Ap 
2po + Ap 



Ap 



u 



4 



+ 4k 2 {Au) 2 p 3 



1,2 



(81) 



(82) 



(83) 



Theorem 6.2 For any integer p there is a solution of the system of equa- 



tions (61) in the sense of the definition 6. Moreover, 

{ p (An) 2 (p + Ap) + ACTAp = 0, 
Act (po + ^Ap) + k 2 Ap = 0, 
v = u + Au + p ^. 



Shock profiles of the considered system ( |65| ) can be found on pictures 
We considered case when p 
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13 and An 



1, 



(Fig. H 

no = 1, k 2 = 0.1. We can also take any real po,o"o (see conditions of the 
theorem |6.2[) but here we took po = 1.1, cto = 0.5 and then we should 
calculate Ap and Act. We consider two different types of shock profiles. 



The first is on the Fig. 21, 021 The second is on the Fig. E3j, 24. 
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Figure 21: Shock profiles, v = 1.1417. Figure 22: Shock profiles, v = 0.35833. 
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Figure 23: Shock profiles, v = 1.1417. Figure 24: Shock profiles, v = 0.35833. 



7 Conclusions and remarks. 

In fact, we considered only special kind of solutions from the sets I and 
J. Moreover, mentioned solutions are "approximate" solutions. It is open 
question about existence of the solution of the Hopf equation in sense of the 
definition 3.3 when p = oo. 

We should notice that there is also a Non- Archimedean approach which 
is developed by V.Vladimirov, I.Volovich, E.Zelenov 



], A.Khrennikov 13]. 
This approach based on p-adic valued distributions and used for the con- 
struction of some models in Mathematical Physics. 

The authors of the papers JjJ ; 0] consider the same equations but they 
speculated a different ideology for generalized solutions and generalized func- 
tions. 

In conclusion we should emphasis that our calculation method looks 
like the Fourier method for linear differential equations but applied to the 
nonlinear equations. Compare the method of mode superposition for a string 
and our method for the shock. Our method allowed to obtain all known 
formulas for the shocks characteristics and, in addition, find a microscopic 
behaviour of shocks in the thin layer with an assumption that the profile of 



25 



the shock can be approximated by the orthogonal system of functions. We 
can use Laguerre functions, harmonic functions or any orthogonal system in 
our calculations instead of Hermite functions. 

Our method one can apply to the problems of hydrodynamics, quantum 
mechanics and non-linear optics. 
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